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Abstract 



The Novello-DeLorenci-Luciane (NDL) field theory of gravitation predicts that 
gravitational waves (GWs) follow geodesies of a modified (effective) geometry with 
a speed lower than the velocity of light. The theory also demonstrates that GWs ex- 
hibit the phenomenon of birefringence, formerly believed to be exclusive of electromag- 
netic waves. Here prospective astrophysical tests of these predictions are proposed. 
I point out that future measurements of gravitational waves in coincidence with a 
non-gravitational process such as a neutrino burst (and likely a burst of gamma-rays) 
may prove useful to discriminate among all the existing theories of gravity. It is also 
stressed that microlensing of gravitational waves emitted by known galactic sources 
(i.e., pulsars) in the bulge, lensed by either the Galaxy's central black hole (Sgr A*) or 
a MACHO object adrift among the Milky Way's stars, may provide a clean test of the 
birefringence phenomenon implied by the NDL gravity theory. 
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1 Einstein's Theory and the Equivalence Principle 

A succesful theory of gravity should be able to correctly predict the way this interaction 
occurs in all process in nature. Einstein's theory of gravitation has till now passed all of 
the tests in this concern. However, it encompasses an implicit statement concerning the 
way gravity-gravity interaction develops when compared to gravity-nongravitational energy 
interactions. General relativity stands on the equivalence principle, which states that any 
sort of matter including massless fields like the photon, interacts with gravitational fields 
fundamentally inasmuch as the same manner. This statement allows to interpret all the grav- 
itational interactions, including gravity-gravity as well (this one having no any experimental 
or observational foundation), as due to changes in the space-time geometry induced by the 
presence of matter fields: g^ u = 7^ + ip^. However, if one dismisses the assumption that 
the gravitational energy should encompass the hypothesis of universality of the equivalence 
principle, i. e., Einstein equivalence principle does not apply to free falling "gravitons", a 
field theory of gravity in which the gravity-gravity interaction occurs in a rather different 
way compared to gravity-nongravity can be formulated JJ]. 

The Novello-DeLorenci-Luciane (NDL) theory of gravitation has been recently introduced [p]]. 
It was shown it incorporates essentially all the ingredients general relativity endowes|l[, and 
in this vein it resembles Einstein theory as far as the first post-Newtonian approximation 
for solar system tests is concerned, and also for the radiative solution up to the quadrupole 
formula level. It has been demonstrated that the most striking prediction of the NDL theory 
is related to the velocity of propagation of gravitational perturbations In Ref. was 
shown that gravitational waves (GWs) travels in the null cone of an effective geometry with 
a speed lower than the velocity of light, the one for GWs to travel in Einstein's theory. 

2 The Novello-DeLorenci-Luciane (NDL) Field Theory 
of Gravity 

In a previous paper M a modification of the standard Feynman-Deser approach of field 
theoretical derivation of Einstein's general relativity, which led to a competitive gravitational 
theory, was presented. The main lines of such NDL approach can be summarized as follows: 

• Gravity is described by a symmetric second rank tensor ip^ that satisfies a non-linear 
equation of motion; 

• Matter couples to gravity in an universal way. In this interaction, the gravitational 
field appears only in the combination 7^ + ip^, inducing us to define a quantity 
9ii,v = Ipu + P'iiw This tensor g^ v acts as an effective metric tensor of the spacetime as 
seen by matter or energy of any form except gravitational energy; 

• The self-interaction of the gravitational field breaks the above universal modification 
of the spacetime geometry. 

*A new more stringent test of the NDL theory predictions concerning the birefringence of the GWs will 
be addressed in section 6||. It is shown there that birefringence of GWs is a peculiar characteristic of almost 
all non- linear theories of gravity, except general relativity, and in particular of the NDL. 
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2.1 Notation and Definitions. 

We define a three-index tensor F a pn, which we will call the gravitational field, in terms of 
the symmetric standard variable ip^ (which will be treated as the potential) to describe a 
spin-two field, by the expression^ F a p^ = |(^[ a;i g] + F^j^), where F a represents the trace: 
F a F a ^ ll/ / y^ L (p,a ^Pcni\v^i^ • 

From the above definition it follows that Fap^ is anti-symmetric in the first pair of indices 
and obeys the cyclic identity, that is: F a[W + F mv = and F afll/ + F^ ua + F vafl = 0. 

The most general non-linear theory must be a function of the invariants one can construct 
with the field. There are three of them which we represent by M,N and W, that is: M = 

* a/3X 

F aixv F a ^, N = F^F^,W = F a(3X F = ±F Q/3A F fluX rf*^. 

We will deal here only with the two invariants U = M — N and W. The reason for this 
rests on the linear limit. Indeed, in order to obtain the standard Fierz linear theory — as 
it is required of any candidate to represents the dynamics of spin-two — the invariants M 
and N should appear only in the combination U. This is the case, for instance in Einstein 
General Relativity theory. 

Under this condition, the general form of the Lagrangian density is given by: L = 
L(U,W), with the gravitational action expressed as: S — J d A x^J— 7 L, where 7 is the 
determinant of the flat spacetime metric 7^ written in an arbitrary coordinate system. 

From the Hamilton principle we find the following equation of motion in the absence of 
material sources: 



L F M^) + l w F 



X(p,u)' 



Lx represents the derivative of the Lagrangian with respect to the invariant X, which 
may be U or W. 



3 Velocity of Gravitational Waves 

The GWs dispersion relation in the NDL theory reads: k fl k u [~f lll/ + A^] = 0, where A Miy = 
2-£^[F^ /3 F ( c t/3 - ) — F^F"], with Lu and Ljju corresponding, respectively, to the first and 
second derivative of the Lagrangian of the theory respect to the invariant U, defined below. 

Thence the discontinuities of the gravitational fields propagate in a modified geometry 
which changes the background geometry 7^ (the Minkowski metric) into an effective one 

C ff = 7"" + A"", (2) 

which has dependence upon the field F^^ and its dynamics. The overall characteristic 
of the new geometry is determined by the non-linear character of the lagrangian on which 
the theory is based. Then the GWs velocity (for a massless graviton) in the NDL reads 

2 We are using the anti-symmetrization symbol [x, y] = xy — yx and the symmetrization symbol (x, y) = 
xy + yx. Note that indices are raised and lowered by the Minskowski background metric 7 M „. The covariant 
derivative is denoted by a semicomma ';' and it is constructed with this metric. 



with the velocity of light c = 1 in geometric units. Here we define Z^ v = F^^FV^s — 
F^F V . In the expression for the velocity of the GWs, Eq.@, the Born-Infeld type Lagrangian 
density 



with b a constant and k = w t N , is the most general functional of the invariant of the 
theory U. The quantity U, the dynamical parameter of the NDL theory, is defined in terms 
of the two fundamental invariants of the theory: M and N. Note that in the linear regime 
C(U) = U . We then obtain the standard weak-field limit as it should be for any massless 
spin-2 theory of gravity, including general relativity. The reader can see Ref.|§ for a more 
detailed discussion of the NDL gravitation. 

Thus a crucial test of the NDL theory, and consequently a potential discriminator among 
the existing theories of gravity, could be an exact determination of the velocity of propagation 
of the GWs themselves. This is an issue which is expected to be accomplished with the advent 
of the new generation of GW detectors such as the interferometers LIGO, VIRGO, GEO- 
600, and the TIGAs resonant-mass omni-directional observatories 0]. Below we suggest a 
prospective astrophysical experimental test of the NDL theory involving the detection of 
GWs in coincidence with a neutrino burst from a supernova explosion, including collapsars 
or hypernovae events. 

We stress that the future detection of the GWs themselves (at least for one detector) 
is unable to provide the looked for discriminator criteria to settle this issue in the light of 
Einstein's and NDL theories. Therefore, a non-gravitational astrophysical or cosmological 
process is called for, and the expected neutrino bursts from both the deleptonization process 
in the supernova core and the gamma-ray burst surge accompanying the GWs in a hypernova 
event may prove useful. 

4 Neutrino-Driven Supernovae and Gravitational- Waves 
4.1 Core-Collapse Neutrino-Driven Explosions 

During the precedent three decades most researchers in supernovae physics have explained 
type-II events as a consequence of neutrinos carrying the huge binding energy of the newly 
born neutron star. Then neutrinos deposit a portion of their energy in a low density region 
surrounding the star's core and a fireball of pairs and radiation finally explodes the remainings 
of the star. In these lines, core-collapse supernovae explosions are one of the most powerful 
sources of neutrinos i/ e , i/^,i/ t and its antiparticles, and likely the sterile one v s . Different 
theoretical and numerical models of type II supernovae explosions || || [7j have estimated 
that 




(4) 
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„ /10 km\ ( M NS \ 

AEtotoi = 5.2 x 10 53 erg ^- 5 

V R ns J \IAM Q J v ; 

are carried away by neutrinos. Almost ~ 10 58 neutrinos of mean energies (10 — 25) MeV 
are released over a time scale of seconds through the process 7 + 7 — > e + + — > v + v. 
Investigations have shown that nearly 99% of the total gravitational binding energy of the 
protoneutron star can directly be carried away by these neutrinos on their diffusion timescale 
At nu ~ 12 s after the core bounce AtcB ~ 1 m s @) §[• The remaining energy being radiated 
in electromagnetic and gravitational waves. 



4.2 Gravitational- Wave Characteristics from Local Supernovae 

On the other hand, during the core-collapse of supernova the time-varying anisotropic dis- 
tribution of density gradients in the proto-neutron star translates into the equivalent of a 
changing quadrupole mass-tensor whose dynamics induces emission of gravitational wave 
bursts [ll|. Because the NDL theory agrees with general relativity upto the first post- 
Newtonian order, we can compute the amplitude of the GW signal as 

= 2GfQ^_ ^ 10 ({-18} { -19}) (6) 

3 c 4 D dt 2 v ; 

for distances as far as the Large Magellanic Cloud D ~ 55 kpc. Here Qij defines the 
mass quadrupole tensor. This signal can last for hundred of milliseconds, with maximum GW 
frequency ~ 1kHz. Since the GWs do not couple to any other form of energy they stream 
away from the SN core whereas ordinary neutrinos in principle do not. This interaction 
induces a time-delay in the neutrino propagation respect to light, or equivalently to GWs in 
the Einstein theory of gravitation. We suggest that such time lag can be used also to test 
the prediction of the NDL theory that GWs travel at a speed lower than the corresponding 
one for light. 



4.3 Neutrinos from GRBs 

Current models of GRBs predict both ultra high, very high|| and high energy neutrinos || 
and ultra high energy cosmic rays emissions which may account for the extra-galactic high 
energy proton flux observed. Next we discuss how the most energetic neutrinos (expected to 
accompany the GWs burst from a collapsar) are emitted according to the GRBs standard 
fireball model. The reader can see Ref.|| for a more complete review of this mechanism. 
In the GRBs fireball picture the detected 7-rays are produced via synchroton radiation 
of ultrarelativistic electrons boosted by internal shocks of an expanding relativistic blast 
wave (wind) of electron-positron pairs, some baryons and a huge number of photons. The 
typical synchroton frequency is constrained by the characteristic energy of the accelerated 
electrons and also by the intensity of magnetic field in the emitting region. Since the electron 
synchroton cooling time is short compared to the wind expansion time, electrons lose their 
energy radiatively. The standard energy of the observed synchroton photons is given by 
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where L 7) 5i defines the energy released in GRBs with L 7 = 10 51 L 7j51 ergs -1 the standard 
luminosity of BATSE observed GRBs, At = lAt ms ms is the typical timescale of variability, 
T = 300r 30 o the Lorentz expansion factor, and £b corresponds to the fraction of energy 
carried by the magnetic field: 47Tr^cT 2 B 2 = 87t£bL, being L the total wind luminosity, and 
£ e the one electrons carry away. No theory is available to provide specific values for both £g 
and £ e . However, for values near the equipartition the model photons' break energy E h v is in 
agreement with the observed one for T ~ 300 and At = 1 ms, as discussed below. 

In the acceleration region protons (the fireball baryon load) are also expected to be 
shocked. Then their photo-meson interaction with observed burst photons should produce 
a surge of neutrinos almost simultaneously with the GRBs via the decay 7r + «-> fi + + <-> 
e + + v e + + z/ M . The neutrino spectrum in the fireball driven explosion follows the observed 
7-rays spectrum, which approximates the broken power- law: cx E@, with (3 ~ 1 for 

low energies and (3 ~ 2 for high energies compared to the observed break energy E? ~ 1 
MeV, where (3 changes. The interaction of protons accelerated to a power-law distribution: 
oc E~ 2 , with the fireball photons results in a broken power- law neutrino spectrum: 

oc E' 13 , with (3 = 1 for E v < E h v , and f3 = 2 for > E b v . Thus the neutrino break 
energy E h v is fixed by the threshold energy of photons for photo-production interacting with 
the dominant ~ 1 MeV fireball photons, and reads 

Elc 5 x 10 14 r 2 30 o eV. 



lMeV 



5 Collapsars, Neutrino and Gravitational- Wave Bursts: 
A Test for the NDL Velocity of GWs 

The just described picture for driving supernovae explosions is by now being considered 
unable to explain the observational fact that some supernovae appear to require more energy 
(an order of magnitude higher) than is provided by the current mechanism based on neutrino 
transport ]ll|. Moreover, the trend in gamma-ray burst (GRBs) modelers is converging on 



a scenario in which a massive presupernova star (and its final explosion as a "hypernova" 



is the leading candidate[|Tl|]. This new paradigm the collapsars: supernovae explosions in 
which a stellar mass black hole, formed previously to the star final disruption, is the central 
engine for the GRBs. This model is supported by the fact that some supernovae have 
been found to be associated with GRBs events. The abrupt fallback (AT acc < 10~ 3 s)Q of 
a surrounding accretion disk, remnant of the failed supernova previous stage, triggers the 
emission of strong GRBs most likely accompanied by GWs and neutrino bursts. In our view, 
this model comprises the necessary non-gravitational astrophysical processes (7 + 1/ bursts) 
through which we can stringently test the NDL theory concerning the velocity of propagation 



3 This timescale will define also the main characteristic frequency of the GW signal emitted. 
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of GWs. For more details on the collapsar mechanism we address the reader to Ref . [TT| , and 
references therein. 

Thus, let us assume for a while that the gravitational radiation (including the GW 
burst produced during the neutrino outburst (TO]) travels at the speed of light. Because of 



the observational evidence that neutrinos actually oscillate [|l3j, which implies they endow 



a mass; and consequently cannot travel at the speed of light, we can use the analogous 
expression for computing the neutrino time delay compared to photons emanating from the 
heavy neutrino radiative decay channel, to estimate their proper time delay with respect 
to the gravitational radiation surge generated at core bounce. Then the time delay for the 
neutrinos (emitted simultaneously with the burst of GWs) to arrive to the neutrino telescope 
is expressed as [OL 0, [UJ 




AT GWW . = 0.515 s — — — , (9) 



where E v represents the neutrino energy, D the source distance to Earth and m v the neu- 
trino mass. Since there is a network (SNEWSQ) of neutrino detectors currently running that 
are sensitive to the prompt core-collapse supernova neutrino bursts in our galaxy JT7[|, which 
can include futurely the new generation of GWs observatories already near completion [|j], 
the appropriate timing of both signals [y + GWs) will provide the time-of-ffight lag in be- 



tween, i. e., the neutrino time delay will directly be stablished by both the observations [pTq] , 
provided the source pinpointing by both capabilities be settled. 

Thus, for a 10 kpc distance, e. g., to the galactic center; for instance, the expected 

neutrino time lag should be: ATgws< >v s = 0.515s, for a (z/ e ) neutrino mass < lOeV, and 

energy < lOMeV, as in SN1987A. Thus the comparison between measured and theoretical 
time-of-flight delay will lead to a highly accurate estimate of the GWs velocity. An inferred 
mismatch between both timescales (expected and measured) may signal that the GWs speed 
as predicted by Einstein theory is not the correct one. This fact would positively point 
towards the NDL prediction as a more plausible explanation, since alternative theories as 
scalar-tensor gravity or other bi-metric gravitational theories predict that GWs travel at the 
speed of light, too. 



6 Birefringence of Gravitational Waves: The Method 
of the Effective Geometry 

Just for later comparison we resume briefly the propagation of electromagnetic waves in a 
non-linear regime. As it was shown |20| the non-linear photon propagates in a modified 
Riemannian geometry characterized by an effective metric g e JJ which is not the background 
one| 
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L Fl ^ - 4 (L FF + VtLpo) F^ X F + (L FG + VLLqq) F\F 



1* Xu 



(10) 



4 The SuperNova Early Warning System. 

5 Although the use of this formulation is not mandatory, it simplifies greatly the analysis of the properties 
of the wave propagation. Besides, we can describe the non-linear photon propagation in a frame in which 
the electromagnetic forces on the photon are eliminated. 
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Here the scalar Q obeys the equation: Q 2 Qi + Q Q 2 + ^3 — 0. The proof of this was 
presented in . 



6.1 The Case of Gravitational Waves: One-parameter Lagrangians 

Our main purpose in this section is to investigate the effects of nonlinearities in the equation 
of evolution of gravitational waves. We will restrict the analysis in this section to the simple 
class of LagrangiansQ defined by L = L(U). 

From the least action principle we obtain the field equation 



= 0. (11) 



~LuF x ^ 

Using the antisymmteric and cyclic properties of the tensor F a> pj (and its covariant 



derivatives) for the discontinuity of the field equation (JllD through the Hadamard's surface 
£ we obtain 

U(af0 ^ + 2 -p— £ Fn{aj3) F = (12) 

where f is defined by f = F a ^ f a ^ - F^ f^. The consequence of such discontinuity in 
the identity (the analogous of the electromagnetic cyclic condition F^-x} = 0) 

F a p \ x + Fp\ v . a + F^ = \{8lW m + 5pW [Xa] + KWwh (13) 

with W a p = F a £ p. £ — F a> p, W\ a /3] = W a p — Wp a , yields after some algebraic steps 
£77^ fc" k v - 2 F afJ/U fp^l k a k 13 + F a ? x f a hp k x + F a fp k a k p = 0. From these equations we 
obtain the propagation equation for the field discontinuities 

LutTKK + ALuu{F^ F\p - F a p u - ^F v )k^ v = 0. (14) 



Expression ([14]) suggests that one can interpret the self-interaction of the background 
field F^ ua , in what concerns the propagation of the discontinuities, as if it had induced a 
modification on the spacetime metric 77^, leading to the effective geometry 

g% = L v rT + 4 L VU (F^ F\p - F aP v - ^ W). (15) 

A simple inspection of this equation shows that in the particular case of the linear the- 
ory the discontinuity of the gravitational field propagates along null paths in the Minkowski 
background. The more general case with two-parameter will be discussed elsewhere [||]. Thus 
the last equation confirms that the propagation of gravitational waves also exhibit the bire- 
fringence phenomenon. 



6 The NDL theory is contained in this class. 
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7 Astrophysical Test for Gravitational Waves Birefrin- 
gence 

The MACHO Collaboration has announced recently that astronomical observations of starfields 
in our galaxy, using the Hubble Space Telescope (HST) and ground based telescopes, have 
provided compelling evidence for the existence of stellar-mass (~ 6 M ) black holes (BHs) 
adrift among the stars comprising the Milky Way||21||. The two observations (1996 and 1998) 
revealed a subtle brightning of a background star produced by the microlensing gravitational 
enhancement of the light it emits due to the passage of an invisible object in between the 
star and Earth. A detailed analysis of the data ruled out white dwarfs or neutron stars 
as the lensing invisible source, and strongly points towards dark stellar-mass objects (i. e., 
black holes) as the magnification sources since ordinary (massive) stars would be so bright 
to outshine the background star. These observations could have been supplemented by the 
discovery of multiple images of the lensed star but unfortunately the HST angular resolution 
is about two orders of magnitude larger than the minimum required for resolving (observing 
the separation of) a pair of images from it induced by the BH bending angle p2| 



a(r 



^ 2 dr , , 
TJ2- 71 - ( 16 ) 



b2 e r 



Here tq is the passage distance from the lensing object, fi and v are the metric fields of 
a static spherically symmetric object, and b is the impact parameter 

b = r e- u{ro)/2 . (17) 

In addition to these effects the starlight (radio waves, for instance) should undergo a 
time-delay respect to a pulse traveling in a region free of gravitation which may be measured 
by precise timing or throughout polarization patterns from the star. This effect, the Shapiro- 
Delay, is due to the light travel through changing gravitational fields. It was also predicted 



to exist for the case of binary radio pulsars. p3|, |24 



In the lines of this microlensing effect of starlight by a BH, analogously a gravitational 
wave (GW) signal from a galactic background source, a compact binary pulsar like PSR 
J1141-6545 (5 hours period), PSR 1534+12 or PSR 1913+16 should also be lensed (splitted) 
when passing near a massive compact dark object such as the MACHO Collaboration BHs. 



Since both theoretical accurate estimates ||25|| and observational statistical inferences [pq] of 
the abundance of galactic neutron star-neutron star binaries and coalescence rates of them 
are more promissing than earlier calculations, the following astronomical configuration looks 
a target to search for. Let us think for a while that a galactic but distant binary radio 
pulsar is aligned with the lensing object (a Schwarszchild BH) and the Earth. A GW pulse 
is emitted from the binary, passes by the lens and is detected at Earth. Then according 
to general relativity both the polarization modes h + and h x of the (linearized) GW signal 
will undergo deflection and time delay when flying-by the lens as in the case for light waves, 
but both will arrive to the detector at the same time, that is, there will be no time lag 
because in GR GWs travel at the speed of light and there is no birefringence effects on their 
propagation. Nevertheless, a phase lag for them in GR is predicted to be exactly 7r/4 radians. 
This dephasing is expected to be measured by the new GWs detectors |27j]. As expected the 
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signal power should be enhanced (enlargening of the GW amplitude) in a forseeable manner 
(For a more extensive discussion of this issue the reader is addressed to De Paolis, Ingrosso 
and Nucita (2000) @). 

Notwithstanding, in the NDL theory of gravity the existence of birefringence of the 
gravitational waves as described above will induce not only a rather different time delay but 
a phase lag too in the arriving GW signals, due to the different velocity of propagation Vk for 
different spatial directions, as showed earlier (see Ref. |T| for further details and definitions). 
This property may be tested with data collected with the forthcoming generation of GWs 
observatories such as LIGO, VIRGO, GEO-600, TIGAs, etc. cross-correlated with data 
from neutrinos, gamma-ray bursts and cosmic rays detector s [J29| . We have shown above 
that each polarization mode of the GW in the NDL theory is velocity-dependent (upon 
direction and magnitude). Then, the radiation component traveling at the lens left-hand side 
(from our point of view) will be affected in a different way compared with the right-handed 
component due to this global birefringence dependence. Thus the detected signals will be 
accordingly time lagged and phase-modulated in a way not mimicking GR, and such effects 
may be measured futurely. The above astrophysical scenario also works for a gravitational 
radiation source at the other side of our galaxy intervened by the Milky Way central black 
hole candidate Sagitarius A*. Moreover, if the lens BH is a Kerr type one then the frame 
dragging (Lense-Thirring effect) induced by the BH spin would dramatically accentuate these 
effects, and it turns out that its observational verification will be a reachable endeavour in 
the days to come. 



8 CONCLUSIONS 

Above we have shown how the almost simultaneous emission of GWs, GRBs and us in a 
single astrophysical event may provide the non-gravitational processes that may turn the 
discrimination between general relativity and the NDL theory of gravity a reachable task in 
the near future. Prospective timing (detection) of such bursts from a unique source on the 
sky may prove powerful to settle the discrepancy between both theories in what concerns to 
the velocity of propagation of GWs. In this sense, the new generation of gravitational- wave 
observatories such as LIGO, VIRGO and GEO-600, together with the SNEWS neutrino 
network and the GRBs new detectors, and potentially the ultra high energy cosmic rays 
observatory AUGER, may prove useful. Moreover, because the neutrino energy can be 
measured by the time it gets the neutrino telescope and the source distance can be reliably 
estimated as discussed above, then from Eq.(|9|) the mass of the neutrino responsible for 
the observed event will be determined or stringently constrained by means not explored 
earlier. This will yield an innovative manner to check the threshold set to the neutrino mass 
by SuperKamiokande neutrino detector contained events. The point here is that despite 
the occurrence of several uncertainties (approximations, etc,) in the derivation of Eq.@, 
the actual detection of both signals by the respectives observatories will render the task of 
constraining the velocitiy of propagation of the GWs and the mass of the neutrino involved in 
the process a feasible one. This is a paramount and inedit manner of weighting the neutrinos 
and measuring the speed of the gravitational waves. 
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